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Abstract 

In this paper, we establish continuity of the path delay operator for networks whose 
arc flows are described by the link delay model presented in Friesz et al. (1993). 

1 Introduction 



One way to model path delay in dynamic networks is that proposed by Friesz et al. ( 19931 ). who 
employ arc exit time functions together with a volume-dependent link traversal time function. 
Such a perspective on path delay has been re ferred to both as t he link-delay model (LDM) 
and as the point- queue model (PQM). Since the lFriesz et al. (|l993T ) paper, conjectures but few 
results about the qualitative properties of the LDM/PQM delay operator have been published. 
One result that is needed for the study of dynamic user equilibrium (DUE) existence as well 
as for analyses of DUE algorithms, when network loading is based on the LDM/PQM, is 
continuity of the path delay operator. 



2 Background and preliminaries 

This paper is concerned with one type of dynamic traffic assignment (DTA) known as continuous- 
time simultaneous route-and-departure-choice dynamic user equilibrium for which unit travel 
cost, including early and late arrival penalties, is identical for those route and departure time 
choices selec ted by travelers be tween a given origindestination pair. Such a mode l is or i ginall y 
presented inlFriesz et al. I (Il993l ) and discussed subsequently bv lFriesz et all d200ll . l20lll201ll ): 
Friesz and Mookheried ()2006l ). 



2.1 Dynamic user equilibrium and delay operator 

We begin by considering a planning horizon [to, tf] C R+. Let V be the set of paths employed 
by road users. The most crucial ingredient of the DUE model is the path delay operator. Such 
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an operator, denoted by 

D p (t, h) for all p G V, 

maps a given vector of departure rates h to the collection of travel times. Each travel time is 
associated with a particular choice of route p G V and departure time t G [to, £/]• The path 
delay operators usually do not take on any closed form, instead they can only be evaluated 
numerically through the dynamic network loading procedure. On top of the path delay oper- 
ator we introduce the effective path delay operator which generalizes the notion of travel cost 
to include early or late arrival penalties. In this paper we consider the effective path delay 
operators of the following form. 

%(t,h) = D p (t,h) + F[t + D p (t,h) -T A ] forallpGP (2.1) 

where Ta is the target arrival time. In our formulation the target time Ta is allowed to 
depend on the user classes. We introduce the fixed trip matrix (Qij '■ (i, j) G W), where each 
Qij G M+ is the fixed travel demand between origin-destination pair G W. Note that Qij 
represents traffic volume, not flow. Finally we let Vij C V to be the set of paths connecting 
origin-destination pair G W. 

As mentioned earlier h is the vector of path flows h = {h p : p G V}. We stipulate that 
each path flow is square integrable, that is 



he{£ 2 + [t , t f }) 
The set of feasible path flows is defined as 



\T\ 



A=\h>0:Yl [ f h p (t)dt = Q ij for all G W \ Q {L\ [t ,t f ]) ]ri (2.2) 

Let us also define the essential infimum of effective travel delays 

= essinf [^ p (t, h) : p G Vij] for all G W 



The f ollowing definition of dynamic user equilibrium was first articulated by iFriesz et al- 
ii 19931 ): 



Definition 2.1. (Dynamic user equilibrium). A vector of departure rates (path flows) 
h* G A is a dynamic user equilibrium if 

h; (t) > o,p g Pij => * p [t, h* (t)] = 

We denote this equilibrium by DUE (^ , A, [to,tf]). 



Using measure theoretic arguments, IFriesz et al.l ( 19931 ) established that a dynamic user 



equilibrium is equivalent to the following variational inequality under suitable regularity con- 
ditions: 

find h* G Ao such that 

E f tf %(t,h*)(h p -h* p )dt>0 }VI(%A,[to,t f ]) (2.3) 
per J t 

for all h G A 
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2.2 The link delay model (iFriesz et all 1 19931 ) 

We shall consider a general network (A, V) where A and V denote the set of arcs and the 
set of nodes, respectively. Additionally, we shall take the link traversal time to be a linear 
function of the arc volume at the time of entry. We describe arc volume as the sum of volumes 
associated with individual paths using this arc, that is 



x a (t) = ^^<5 a pX^(i) for all a € A 

pev 



(2.4) 



where Xa(t) denotes the volume on arc a S A associated with path p 6 V, and V is set of all 
paths considered. Each path p € V is taken to be the set of consecutive arcs that consitute 
the path; that is 

p = {ai, . . . , a m ( p )} 
Furthermore we shall make use of the arc-path incidence matrix 



A = (8, 



ap J 



where 



'ap 







if arc a belongs to path p 
otherwise 



We also let the time to traverse arc for drivers who arrive at its entrance at time t be 
denoted by D ai (x ai ). 

It will be convenient to introduce cumulative exit flows V ai (-), Vj^(-): 



Jto Jto 



ds, 



p £V, cii £p 



where the notation employed has obvious and conformal definitions relative to that introduced 
previously. The following diff erential alge b raic e quation (DAE) system is another version of 
the DAE system employed by lFriesz et al.1 (120111 ): 



Vl^t) = Vl(t + D ai {X ai (t)))i 
Furthermore, we make the following definitions 



vUt) = h p (t), 



1, . . . ,m(p) 



h p (s) ds 



(2.5) 
(2.6) 



(2.7) 



to 



where h p , also known as the path flow, is the rate of departure from the origin of p £ V . It is 
also conventional to introduce the link exit time function 

T a (t)=t + D a [x a {t)} 

for each a € A. 

The next theorem, from lFriesz et al.l (| 19931 ). presents an important property of linear link 
delay functions: 

Theorem 2.2. For any linear arc delay function of the form D{x) = ax + (3, a, (3 > 0, the 
resulting arc exit time function r is continuous and strictly increasing and hence the inverse 
function r" 1 exists. 

Proof. See Theorem 1 of lFriesz et al. (|l993h . □ 



3 



3 The main result 



The following is a statement of our main result: 

Theorem 3.1. Consider a general network (A, V), where arc dynamics are governed by the 
link delay model, assume the link delay function for each a G A is afflne. That is 

D a [x a {t)\ = a a X a (t)+p a 

V\ 



where a a G 'St] r and (3 a G 8^4.4.- Then the effective delay operator from A into \^C 2 ([to, tf] 
h G A — > h) is a continuous map. 

Remark 3.2. In Zhu and Marcotte hood ), the authors show that the effective delay operat 



or 

is weakly continuous when the LDM is employed, under the restrictive assumption that the 
path flows are a priori bounded from above. Such assumption is dropped in our result; we also 
assert strong, not weak continuity. 

The above theorem tells us there are no jump or other kinds of discontinuities of the path 
delay operator. Such an analytical result is crucial for the study of dynamic user equilibrium 
existence as well as for analyses of DUE algorithms, when network loading is based on the 
LDM/PQM. 

3.1 Proof of the main result 

Now we present the proof of Theorem 13.11 

Proof. We begin by showing that given a converging sequence hP^ in the space A C (^£+([to, t 
such that 

||/i (n) - h\\ c2 — ► n — > 00, (3.8) 

the corresponding delay function -D p (-, hS n '\ converges uniformly to /i) for all p G V. 

This will be proved in several steps. 

Part 1. First, let us consider just one single arc, and hence omit the subscript a for brevity. 
Assume a sequence of entering flows {u^} n >\ converging to u in the C 2 ([to, tf]) space; that 

\\u {n) -u\\ 2 = Qf' (u^(t)-u(t)ydtj — > n —> 00 (3.9) 
Dehne the cumulative entering vehicle counts 

Wi) ± /%(->(.,* „ > 1 

t 
to 

Then we assert the uniform convergence UW — ► U on [t , t f }. To see this, fix any e > 0, in 
view of (I3.9p . choose > such that for all n > N 

\\u^ - u\\ 2 < e 



to 



t G [t , t f ] 
U(t) = I u{s)ds 



\v\ 
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According to the embedding of /^([to, tf]) into C 2 ([t , tf]), we deduce for any t £ [to, tf] that 



\U (n \t) - U(t)\ 



u^\s)ds- / u(s)ds 

to J to 



< \\vF> - < (t -t f ) 1/2 \\u {n) -u\ 



< (to-tf) 1 ^ 

The preceding shows the uniform convergence — on \tn t f}. 

Part 2. We adapt the recursive technique devised in iFriesz et all (Il993h . Let X^(-), n > 1, 
and X(-) denote the arc volumes corresponding to U^ n \-), n > 1 and U(-), respectively. 
Assume, without loss of generality, that 



X in \t c 



0, 



X(to 







and that, for the flow profile U(-), the first vehicle enters the arc of interest at time to- I n 
addition, let t\ denote the time that first vehicle exits the arc of interest. By definition 



t x = D(0) = p 

For all t £ [to, ti], since no vehicle can exit the arc before time ti, we have 

X (n) {t) = ^(n) (t) X{t) = U[t) tG[t0jtl ] 



(3.10) 



(3.11) 



For each flow profile U^ n \ n > 1, denote the exit time function restricted to [to, ti] by r[ n \-); 
under flow profile U, denote the exit time function restricted to [to, ti] by ti(-). Then 



r W(t) = t + D(xW(t)) = t + aU( n \t) + p, t€[to,h] 
n(t) = t + D(X(t)) = t + aU(t)+/3, te[t ,ti\ 



(3.12) 
(3.13) 



We conclude that rj; 



(n) 



T\ uniformly on [to, t\\. Now let 
t 2 = mfr^Cti) < r x (t x ) 



(n)\-l 



, n > 1, 



Fix 5 small enough and call t 2 = t 2 — S. See Figure [TJ By Theorem 12. 2[ (r-[ 

and r-j -1 are well-defined, continuous and strictly increasing. We claim that {(t^) 1 } n >i 
uniformly converges to Tf 1 on [t%, t 2 ]. To see this, we need to extend the arrival time function 
T\ and to the interval (— oo, to)- Because no vehicle is present during (— oo, to), it is 
natural to assign 

n(t) =t + p, rW(t) = t + p 

This means, if an infinitesimal flow particle enters the arc at t 6 (— oo, to), its travel delay 
will always be (3. Fix any e < t\ — r ] ~ 1 (t 2 ), and consider the following quantities: 



A: 



inf 



l n (t)-n(t-e)} A+ = inf l n (t + e)-n(t)} (3.14) 
L J teltotT- 1 ^ L J 
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Since the infimum of a continuous function on a compact interval must be obtained at some 
point t G [to, r-j - (£2)]} we conclude A~, A+ > by the strict monotonicity of t\ established 
in Theorem E2J 



According to the uniform convergence t\ 
such that as soon as n > iV, we have 



(«) 



ri on (—00, ti], there exists some iV > 



4 n) (t) < n(t) + A~/2 r[ n> (t) > n(t) - A+/2 for all t G [-00, tj 



(3.15) 



For any s G [ii, £2], we have r x G [to, T\ (^2)]- Therefor, for all n > N, in view of (I3.15|) 
and (|3.14p . we have 



e) < n^^-^+Aj^ < r^rf 1 ^)) -A-/2 = S - A"/2 



r[ n \r 1 - 1 (s)+e) > n (rf 1 (a) + e) - A+/2 > n (rf 1 (,)) + A+/2 = » + A+/2 

By the Intermediate Value Theorem, there exists some t* G [rf ^s) - e, r" 1 ^) + e] with 
t[ (t*) = s. In other words, we know 

\{ri n) r\s)-rrHs)\ = [ 



t x \s)\ < e, 



for all n > N 



This finishes our claim that (ti ) 1 — > r x 1 uniformly on [ti, tj\. 

Let T2 (•), n > 1 and T2(-) be the exit time functions for commuters entering during the 
interval [t\, £2]; corresponding to entering flow profiles U^ n \-), n > 1 and J7(-), respectively. 
Then for each t G [ti, ^2] , we may state 



rJ B) (t) 
r 2 (t) 



t + a 
t + a 



+ /3, n > 1 



[^-[/(rf^)) 



(3.16) 
(3.17) 



Now we make the claim that [/( n )^(r^) 1 (t)J — > ^( T i H*)) uniformly on [ii, 43]. Indeed, 
for any e > 0, there exists iVj. > such that, if n > Ni, we have 

|^7 {n) (t) - C/(t)| < e/2, for all i G [i , */] 
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Moreover, the functions (t^) 1 , n > 1, and r-j -1 have a uniformly bounded range on [ti, £2], 
namely [to, ti]. By the Heine-Cantor theorem, {/(•) restricted to [to, ti] is uniformly contin- 
uous, which means we can find 80 > such that, for any si, S2 G [to, ti] with |sj — S2I < <5o, 
the following holds: 

\U( Sl )-U(s 2 )\ < e/2 



By uniform convergence of (rj 
have 



(n)\-l 



T t , we may choose A2 > so that, for n > N2, we 

M n) )-\t)-rr\t)\ < S 
Thus we deduce that, given n > max{iVi, A^}, for any t S [ti, t2], the following is true: 

^ ) (tf ) )' 1 w)-^r 1 w; 



< 



^((rj"')" 1 ^) 



1 



(^l + l^^)" 1 ^))-^ 1 ^)) 



< e/2 + e/2 



^((r^)" 1 ^)) — »• ^(rf 1 ^)) on [t a , t 2 ], and our 



This shows the uniform convergence 

claim is substantiated. As an immediate consequence of (|3,16p and (|3.17p . converges to 
r 2 uniformly on [ti, £2]. 

Part 3. We now proceed by induction as follows. Choose any v > 2, and call 



tu+i — t v+ \ — 5 



where the constant 5 is the same as what was used in Part 2. Using the induction hypothesis 
that t^ 1 converges to t v uniformly on \t u -i, tu], we show the following uniform convergence 
on [t„, t v+1 ]: 

The proof is similar to what has been done in Part 2. Now introduce t^™^-), n > 1, and 
r !/+i(')> which are the exit time functions corresponding to U^(-), n > 1 and U(-) respectively; 
and they are restricted to the time interval [t v , t v +i]. Similar to results fj3. 161) and (13.17p . we 
have for t € [t u , t u +i], that the following holds: 



r&(*) 

Ti/+l(t) 



t + a 
t + a 



C/(«)(t)-C/(«)(( ri ( n) ) -1 (t))l +/3 n > 1 
I7(t) - ^(r-Ht))! +/3 



(3.18) 
(3.19) 



It can be shown as before that [[tv ) X ) — )• C/(r^ x ) uniformly on [t u , t v +i]. Therefore 
r^"\ — > t v+ i uniformly on \t u , t^+i]. This finishes the induction. 

Part 4. We now have obtained a sequence {[t v , tv+i]) u> $ of intervals. On each interval 
[t u , t u +i], the uniform convergence 



An) 
r u+l 



T v +1, 



holds. Notice that, by construction, t y+ \ — t v > (3 — 5 > 0, for all v > 1; therefore the interval 
[to, tf] can be covered by finitely many such intervals. As a consequence, we easily obtain the 
uniform convergence 

of r (n)(.) T (.) on 

the whole of [to, tf], where r^ n \-) corresponds to 

U( n \-), n>l and r(-) corresponds to U(-). 
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Let (•), n > 1, and v(-) be the exit flows of the single arc and then define the cumulative 
exit vehicle count 

V(t) = f v(s)ds, V {n \t) = f t v ^(s)ds, 

J tg J to 

It then follows immediately from the relationships 

V(t) = U^-^t)) V^ n \t) = l/W^rW)" 1 ^)) 
that converges uniformly to V(-) on [to, tf]. 

Part 5. Consider a general network (A, V) with a converging sequence — > h in A C 
(£ 2 ([t , */])) • Define for p £ V the following: 

fl-W(t) = fhf\s),d8, H p (t) = fh p {s)ds 
J to J to 

Then the fl£ n) (-) converge uniformly to H p (-). For each arc a £ A, where the notation 
employed has an obvious meaning, the cumulative entering vehicle count lli n \-) is given by 

P a'eX(o) 

In the above, the first summation is over all paths that use a as the first arc; and, in the 
second summation, 1(a) denotes the set of arcs immediately upstream from arc a. A sim- 
ple mathematical induction with the results established in previous steps yields the uniform 
convergence 

E#°(t) — ► U a (t), y("'(t) — > K(t), rW(t) — ► r a (t), for all a G .4 

where r a (-) is the exit time function of arc a. Thus, for each path p G V, the path delay 
D p (; hW) also converges uniformly to D p (-, h) since it is a finite sum of arc delays. It remains 
to show that the effective delays obey h^) — > h^) uniformly on [to, tf]. Recall 

%(t,h) = D p (t,h)+T(t + D p (t,h)-T A ) 

Notice that J- is uniformly continuous on [to, tf] by the Heine-Cantor theorem; this means, 
for any e > 0, there exists a > such that whenever |si — S2 1 < cr, we have 

\T( Sl ) - T{s 2 )\ < e/2 

Moreover, by uniform convergence, there exits N > such that, for all n > N, we have 

\D p (t, h {n) ) - D p (t, h)\ < min{cr,e/2} for all t G [t , tf] 

We then readily deduce that, given n > N, the following holds for all t G [to, tf]: 

%(t, h,W)-* p (t, h) 

< D p (t, h^)-D p (t, h)\ + \?(t + D p (t, hW)-T A ) -F(t + D p (t, h)-T A ) 

< e/2 + e/2 = e 
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Final Argument. Finally, uniform convergence on the compact interval [to, tt\ implies con- 
vergence in the C? norm: 



£ f ($ p (t, ~%(t, h)y dt — > o, 



n 



oo, 



p G V 



(3.20) 



Summing up (I3.20P over p E V, we obtain the convergence *(-, /i (n) ) — ► *(•, h) in the 

/ \ \~P\ 

Hilbert space [C 2 ([t , t f ])\ . □ 



4 Concluding remarks 

Dynamic traffic assignment differs from static traffic assignment in that path delay does not 
enjoy a closed form. In fact the path delays needed for the study of dynamic user equilibria 
(DUE) are operators that may only be specified numerically. Moreover, such path delay 
operators are based on the specific model of arc delay employed for network loading. We 
have considered in this paper path delay operators for the network loading procedure that 
is endogenous to the Friesz et al. (1993) which has been referred to as the link delay model 
(LDM) and also as the point queue model (PQM). We have shown that LDM/PQM path delay 
operators are continuous under the v ery mild a s sump ti on that the link tr aversal time function 
is affine. Combined with the results, Browder ( 19681 ): Han et al. ( 2012al bl 3). an increasingly 
clear understanding of DUE based on different network performance models is emerging. 
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